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Abstract 

We determine the low-energy spectrum and the eigenstates for a 
two-bosonic mode nonhnear model by applying the Inonii-Wigner con- 
traction method to the Hamiltonian algebra. This model is known to 
well represent a Bose-Einstein condensate rotating in a thin torus en- 
dowed with two angular-momentum modes as well as a condensate in a 
double-well potential characterized by two space modes. We consider 
such a model in the presence of both an attractive and a repulsive bo- 
son interaction and investigate regimes corresponding to different val- 
ues of the inter-mode tunneling parameter. We show that the results 
ensuing from our approach are in many cases extremely satisfactory. 
To this end we compare our results with the ground state obtained 
both numerically and within a standard semiclassical approximation 
based on su(2) coherent states. 

PACS: 03.75.Fd, 03.65.Sq, 03.75Lm 



1 Introduction 

The dynamics of a bosonic fluid rotating within a thin torus and, particu- 
larly, the study of the properties relevant to its weakly-excited states have 
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received recently a large attention PP-jH] due to the rich phenomenology that 
characterizes such a system. For example, the quantization of fluid circula- 
tion is shown [3^ to disappear whenever the physical parameters cause the 
hybridization of condensate ground state over different angular momentum 
(AM) states. A similar effect is found in the mean-field dynamics of the 
condensate wavefunction on a circle |4j, where the circulation loses its quan- 
tized character when the system is in the soliton regime. The rotating fluid 
exhibits low-energy AM quantum states (corresponding to the presence of 
plateaus of quantized circulation) that determine the hybridization effect by 
a suitable tuning of the model interaction parameters jSj- In the simplest 
possible case, the model exhibits two momentum (bosonic) modes associated 
to two AM states (the ground-state and the flrst excited state) of the fluid. 
An almost identical model [01-53 been studied thoroughly in the recent 
years within Bose-Einstein condensates (BEC) physics, where a condensate 
is distributed in two potential wells that exchange bosons via tunneling ef- 
fect. The two-well model H = U {nl + n1) — A{nQ — ni)/2 — Vo{aQaf + aiOQ ), 
where Oq, ai are bosonic space-modes and rii = afai, displays hybridized 
states when the well-depth imbalance vanishes (A = 0). 

For both models the energy regime of interest is that corresponding to 
the ground-state or to weakly excited states. In this respect, many authors 
have tried to develop approximation schemes able to provide a satisfactory 
analytical description of the low energy spectrum and of its states. The non- 
linear character of the model Hamiltonian entails a difficult diagonalization 
process unless one resorts to numerical calculations. In this case the exact 
form of the spectrum is obtained quite easily. However, for N-well systems 
such as condensate arrays described by the Bose-Hubbard model, Josephson- 
junction arrays and, in general, N-mode bosonic systems ^1], ^H], the exact 
diagonalization requires a computational effort rapidly increasing with N. 
This motivates the interest in developing effective, analytical approximation 
methods able to solve the diagonalization problem. 

The present work has been inspired by papers [2] and [3] where, among 
the variuos issues considered, the structure of the ground state of a ring 
condensate (within a two-AM-mode approximation of the bosonic quantum 
field) has been studied. As to the closely related two-well boson model, the 
same problem has been investigated in ^H] within the hermitian phase op- 
erator method. In order to obtain a satisfactory description of the system 
ground state as well as of the weakly-excited states for the two-mode model, 
we implement, in the present paper, an algebraic approach based on the 
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Inonii-Wigner contraction method ^7]. This method allows one to simplify 
the algebraic structure of the Hamiltonian reducing the latter in a form apt 
to perform a completely analytic derivation of its spectrum. A well defined 
limiting procedure, mapping the original Hamiltonian generating algebra to a 
simpler algebra, often succeeds in reducing the nonlinear terms to a tractable 
form. These terms, originated by the boson-boson interaction and thus oc- 
curring in any model inherent in BEC dynamics, are known to make the 
Hamiltonian diagonalization a hard task. Such a technique and the effect of 
simplifying the algebraic structure of model Hamiltonians, has found a wide 
application in many fields of theoretical physics. It is well illustrated, e. g., 
in reference |^ where it is applied to study collective phenomena in nuclear 
models. 

The contraction-method approach (CMA) -namely the contraction pro- 
cedure and the ensuing approximation of weakly excited states- works well 
for the spectrum sectors where the energy levels are close to the minima and 
the maxima of the classical Hamiltonian and thus seems suitable for studying 
the low-energy regime of two-mode nonlinear models. The results obtained 
within the CMA in sections |21 and 01 will be compared both with the exact 
spectrum calculated numerically and with an alternative approch based on 
the coherent-state semiclassical appproximation (CSSA) reviewed in section 

m 

We consider interacting bosons with mass m whose boson-boson in- 
teraction can be either attractive or repulsive. These are confined in a nar- 
row annulus whose thickness 2r is much smaller than the annulus radius 
R. Bosons are also acted by an external potential which causes inter-mode 
tunneling. Particularly, the rotating fluid with attractive interaction can be 
shown to be equivalent to the two-well model of repulsive bosons introduced 
previously. In the coordinate frame of the potential rotating with angular 
velocity u and with z axis parallel to total angular momentum Ltot = Lz, 
the bosonic-field Hamiltonian reads 

r ^ I 
Hbf = J f/'r^J 



uLz + Vext{r] 
Im 
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where ?/'r = V'(r) {i^t) is the destruction (creation) boson field operator 
at r. Vext is the confining potential. At low temperature, the interaction 
between dilute bosons is well represented by the Fermi contact interaction 
which entails the standard approximation U{\v — s|) ~ {ATxTi^a/m)5{\r — s|), 
where a is the s— wave scattering length [2]. 
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1.1 Two-mode approximation 

The two-mode approximation involves only the first two states of AM, with 
eigenvalue equations L^-j/'olr) = and Lzipiir) = htpilr). Field operator tp{r) 
in the two- mode basis of is thus written as V'(r) — aoV'o(r)+«iV^i(i')5 where 
ao, ai are bosonic operators and the validity of the two- mode approximation 
requires the condition < u; < 2uJc (greater angular velocities would involve 
other angular- momentum states). Within such an approximation [HI El and 
considering a thin torus (r << R) H^f reduces to ^ H = gijiQ + n\ — riQ — 
ni + 4?io^i)/2 — A/ini/2 — Vo(a^ao + '3o'o,i), where Ui = afai, A = 2h{uj — Uc), 
while ujc = h/(2mR'^), g = 2^^a/(mi?7rr^), and Vq are the critical angular 
frequency, the mean interaction energy per particle, and the asymmetry of 
potential V^xt = Vo(e*^ + e*^), respectively. In the Schwinger picture [7j of 
algebra su(2) H further simplifies becoming, up to a constant term, 

H = -gJi-2VoJi-AJs, (1) 

where J3 = {ni — no)/2, Ji = (J+ + J_)/2, J2 = (J+ — J_)/2z and J+ = afao , 
J_ = Such generators satisfy the commutators [Jr, Js] = i^rsvJv {^rsv 

is the antisymmetric symbol) and commute with the total boson number 
operator rii + Uq {rii = afai) whose eigenvalue is connected with the 
su(2)-representation index J hj J = 2N. In such a scheme, the AM states 
are defined by 

I J; m) := |no) ® |ni), rii = J + m, tiq = J — m, 

where the Js-basis states satisfy the eigenvalue equations J3I J; m) = m\J;m), 
and J4\J;m) = J\J;m), the index J being the eigenvalue of J4 = (ni -|- 
no)/2. The positive (negative) sign of g in model^implies that the effective 
interaction between bosons is repulsive (attractive). The conditions of weak 
asymmetry and interaction ensuring the validity of the model [2] are given by 
I Vol <^ huc, and \g\ hujc- The simple spin form of Hamiltonian^ evidences 
how the attractive model [g < 0) coincides with a (repulsive) two-site Bose- 
Hubbard Hamiltonian [71 modeling two potential wells of different depth 
that share N = 2J bosons and exchange them via tunnel effect. The A^- 
boson physical states can be written as lip) = J2m=-j ^mlJ^'m) , while the 
Schrodinger equation {ihdt — H)\tp) = can be expressed in components as 

ihXrn = [—g^^ — rnA] Xm — Vq 



m—1 
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once the symbol R'^ = [(J+m)(J — m + has been defined. It worth not- 
ing that the study the algebraic structure characterizing the second-quantized 
Hamiltonian for a condensate trapped in two potential wells has received a 
large attention in the literature. In the seminal work ^3] and in reference 
|15j . in particular, such Hamiltonian has been shown to reduce, within a 
standard mean-field approach, to the sum of mode Hamiltonians describing 
the momentum conservation in the presence of inter- well boson exchange due 
to the tunneling. Each mode Hamiltonian is written in terms of operators 
Ofc, a_fc (±/c are the momentum modes) and can be reformulated as a linear 
combination of su(l,l) generators. In model ^ the momentum conservation 
is explicitly violated since one of the mode takes into account the fiuid ro- 
tation. This fact entails that the previuos Schwinger realization of algebra 
su(2), rather than the algebra su(l,l) connected with the momentum con- 
servation, characterizes the system. 

In our analysis the dimensionless mean-value per boson of the angular 
momentum (l^) = {Lz)/hN (where the notation (A) = {tp\A\tlj) has been 
introduced) represents an important quantity. The angular momentum, in 
fact, expressed as 

relates the macroscopic behavior of the rotating condensate to the minimum- 
energy state properties through the ground-state components X^- In the 
sequel we consider the spectral properties of model [T] both in the attractive 
case {g < 0) 

Ha = \g\ Jl - 2K)Ji - AJ3, (3) 
and in the repulsive case {g > 0) 

Hr = -{\g\Ji + 2VoJi + AJ;). (4) 

It is worth noting that the study the ground-state properties of the re- 
pulsive case is closely related to the study the maximum-energy state for 
the attractive Hamiltonian. In fact, after the substitutions Vq — >■ —Vq and 
A —A, the repulsive Hamiltonian is identical to the attractive one up to 
a factor (—1). Since these two changes can be effected in a unitary way by 
means of transformations e"''*'^"'^ JiC"*'^"'^ = — Ji, and e'^^'^'^^J^e'^'^'^^ = — J3, 
respectively, the spectra of Hr and Ha turn out to satisfy the equation 
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spect[Ha(yo, A)] = —spect[Hr{—Vo, —A)]. Concerning the parameter A of 
Hamiltonian ^ we note that the constraint < u < 2uJc, imphes the in- 
equahty —2hwc < A < 2hwc- The definition of the further parameters 
7 = J\g\/2hwc, T = Vo/J\g\, allows one to better characterize the regimes of 
the rotational dynamics as well as the conditions of validity of the present 
model. Parameter 7 (representing the ratio of the self-interaction energy 
per particle to the single-particle energy-level spacing) should satisfy the in- 
equalities 27 << J, T « 1/27, owing to the conditions \g\ << TiUc and 
Vo << fiujc, respectively. Both these conditions can be satisfied if J = N/2 
is not excessively large Moreover, parameter r = Vo/(J|5'|) allows one to 
distinguish, in both the attractive and repulsive case, three regimes: 

• Fock regime, where \g\ ^ VqJ entails r ^ 1/^/^, 

• Josephson regime, where Vq/J ^ \g\ ^ VqJ entails 1/J^ ^ r ^ 1, 

• Rabi regime, where \g\ <^ Vq/J entails r ^ 1. 

We note that the condition of weak asymmetry | Vo| <^ Tiujc given by r ^ I/27 
appears to be compatible with the first two regimes and with part of the Rabi 
regime. 

2 The Inonii - Wigner contraction in the at- 
tractive case 

We introduce a simple algebraic approach for studying the low-energy spec- 
trum of Hamiltonians 01 and |3] for large J whose essence consists in sim- 
plifying the nonlinearity due to the term J|. The Inonii- Wigner contrac- 
tion ^19j supplies a method for mapping some given algebraic structure in 
a new one, as the result of a singular limiting process. The contraction 
is realized by defining a set of new operators hi as linear combinations 
hi = ail + T^kCikgk of the generators gk of a given algebra (identified by 
its commutators [gr,gs] = £rskgk) and of the identity operator I. Selecting 
an appropriate parametrization Cjfc(x) of the linear- map coefficients, the con- 
traction enacted by means of the limit x ^ is able to generate the new 
algebraic structure [hi,hj] = Cijkhk whose structure constants {cijk} differ 
from the original ones {Srsk}- For the algebra su(2) the contraction of the 
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\9\ Jl - 



algebra mapping is driven by x = 1/ ^/j (with J oo) and generates, in this 
limit, the harmonic oscillator (namely the Heseinberg-Weyl) algebra [20] • 

The classical study of attractive {g < 0) Hamiltonian Ha 
2Vo Jl— AJ3 developed inlXldemonstrates (see formula IH^j) how Ji ~ + J, J 3> 
I J2 1, 1^3 1 — 0, at low energies. This suggests the correct way to implement the 
contraction scheme. In the present attractive case we can build the following 
transformation 



hi = Jl - //x^, /i2 = XJ2, hs = XJ3, 



(5) 



where J4 = J I. The Inonii-Wigner contraction is realized when such a x- 
dependent transformation is considered in the (singular) limit J = l/x^ — *■ 
00. In this case the objects {J,} (with i = 1,2,3,4), defining algebra u{2), 
transform into the new objects {hi, 1} (with i = 1,2, 3) that satisfy the 
following commutation relations: 

[h2, hs] = i{x^hi + /) il, [hi, /i2] = x[Ji, J2] = ih^ , (6) 

[hi, hs] = x[Ji, J3] = -^h2 , [h,,I] = 0. (7) 

In the limit x = — ^ 0, the latter reproduce the commutation relations of 
Weyl-Heisemberg algebra: [q,p] = i, [n,q] = —ip, [n,p] = iq, n = (g^ + p^)/2, 
thereby suggesting the identifications hi = —n, h2 = —p, h-^ = q. By 
combining the latter with definitionsElwe find that the contraction gives Ji — >■ 
J — n, J2 — > —\fjp, J3 — > \fjq. Correspondingly, Hamiltonian becomes 
iia = \g\Jq^ + 2Von - 2VoJ - A^/jq, which, by defining = [1 + l/r]i/^ 
and Q = q — X with x = V^A/2Vof2^, and r = Vq/ J\g\, reduces to the form 



Ha = Vn 



p^ + i]2g2 _ 2 J 



(8) 



Since p^ + Vt^Q^ = 2Q{n + 1/2) is diagonalized by the harmonic-oscillator 
eigenstates ^n{Q) = {Q\En) = A^^e-^^'/^^^ (v^Q) , the eigenvalues of 
Hamiltonian |H1 are found to be 



2VL{n + l/2) -2J 



(9) 



The corresponding eigenvalue equation Ha\En) = En\En) in the J3 ba- 
sis, where \En) = SmX„(m)| J, m), can be written as T,rn{Ha)£m^n{^) = 
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EnXn{i) with {Ha)im = {Jy^lHalJyfn). Ill the hmit J >> 1, equation 
Js\J,m) = m\J,m) is replaced by q\J,m) = {■m/^/j)\J,m). Therefore 
the eigenvalue m/^/j can be seen as a continuous variable which natu- 
rally identifies with the variable q ~ J^/y/j used within the approxima- 
tion scheme just discussed. The component version of the eigenvalue equa- 
tion for Ha then reduces (see reference |^ for details) to the equation 
Ha{QiP)'^n{Q) = En'^niQ) solvcd abovc. Compoueuts Xm{En) thus ap- 
pear to be given by Xm{En) = "^niQ) that entail the explicit expression for 
the eigenstates 

= E^X^(E„)|J;m), X„(E„) = N^H^{VnQ)e~'''^"/^ (10) 

with Q = m/y/j — x- The normalization constants are determined 
through the condition {En\En) = 1 implying that 



J 

E 



^miEn) 



30 /V2 

dq^H, 
oo y J 



Vniq-x) 



(11) 



where ± J has been replaced with ±oo. Such an approximation is acceptable 
until the condition 



2n/VL 
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-evinced from the interval containing the Hermite-polynomial zeros- is ful- 
filled. Excluding the case r ^ 1, this condition is always valid provided 
n « J. Thus constants are given by = [( Jfi)^/V(^^^^2"n!)]^/2^ 

Another important check concerns the possibility of considering m/ \/j as 
a continuous variable. The characteristic scale is established by the gaussian 
deviation 



'2/f2 which must be compared with the smallest variation 1/v J 
of q. The resulting condition l/-\/j < 



2/f2 can be written as 



Vo + J\9\ 



1/2 



2J 



T 



T+1 



1/2 



While in the Rabi and Josephson regimes (1/J^ << t) the latter is fully 
satisfied, in the Fock regime, where r << 1/J^, such condition is violated. 
We notice that for r ~ 1/J^ (namely \g\ ~ JVq) a unique component Xm 
appears to contribute to states \En) since the gaussian amplitude becomes 
very small. For example, in the case of the ground-state one has 



\J,m) ~ Nqc 



(77-^)' 



I J, m*) 



(13) 
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m/J m/J 



Figure 1: In both panels, grey (dark) diamonds describe the ground-state 
components X^, obtained from formula^! for = 20 (A^ = 40) within the 
contraction-method approach (CMA). The edges of the grey/dark piecewise 
hnear curves represent components X^'s calculated numerically. Left panel: 
Josephson regime in the attractive case with r = 0.02, v = 0.8. CMA com- 
ponents X^s and X^'s calculated numerically are almost indistinguishable. 
Right panel: Attractive case with r = 1.0 (transition point from Josephson 
to Rabi regime) and i> = 0.8. No difference is visible between CMA X^s 
and Xm's calculated numerically. 



where m* is the integer closest to ^/Jx — ^/2|5'|- Nevertheless, in the special 
case when A/2\g\ = + 1/2, the two states \ J,m^:) and | J, m* + 1) equally 
contribute to \Eq) which is given by 

l^o) ^ Noe-^{\J,m^) + |J,m, + 1)). (14) 

To summarize, we note how the ground-state |£'o) is essentially formed by 
a unique component corresponding to \J,m^) in the whole parameter range 
— 1/2 < A/2\g\ < + 1/2. The resonance of the system between two 
equivalent states crops up whenever A/\g\ assumes integer values given by 
A/|5f| = 2m + 1 with —J<m< J. Such condition can be implemented by 
varying A with \g\ = const thus leaving Q unchanged. 
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2.1 Comparison of different regimes 

For r > 1/J^ (Rabi and Josephson regimes), one easily calculates the di- 
mensionless mean AM per boson (Z^) based on state \Eq), as given by for- 
mula ^1 and exploiting the normalization integral ^2 Recalling that (J3) = 
J2i=-jmX'^{Eo), one finds 

where (J3) matches exactly formula ESI obtained in the classical study of the 
attractive model. This result cannot be used in the Fock regime where the 
ground state has, at most, either one or two dominating components. In the 
other two regimes, the second of equations El entails the further consistence 
condition 

- 1 < {J3)/J = rA/[2Vo(l + r)] < +1 , (16) 

which has to be verified in each regime. In view of the condition KJs)] << J 
required to implement the contraction procedure, formula El should be im- 
posed in the stronger version |rA/[2Vo(l+T)]| << 1. However, the numerical 
(exact) determination of the ground state for various choices of parameters 
reveals that our approximate procedure works well also in the case when 
|rA/[2Vo(l +t)]| is not particularly small. 

Fock regime. The main feature of this case (r ^ 1/^^) is that the mean 
dimensionless AM per boson is a step function of A (as to this well-known 
effect see, e. g., reference P|). If one simplifies the form of states ITI^andlT^bv 
setting |£'o) = \ J,m) and \Eo) = (|J, m) + |J, m + l))/\/2 in correspondence 
to the appropriate values of A, the dimensionless mean AM per boson is 
found to be 

1 V 1 fTT- 1 

2 2J' ^ ' 2 2J 4J' 

for m — 1/2 < A/2\g\ < m + 1/2 and A/2\g\ =m± 1/2, respectively, corre- 
sponding to the two choices of the Fock ground state \Eo). This illustrates 
the AM step character (related to the Hess-Fairbank effect) as well as its 
"singular" behavior when A/2\g\ = m ± 1/2. Notice that considering the 
simplified form for \Eq) is equivalent to assume the net predominance of one 
or two components. The results just found are consistent with the limit 
r ^ 0, where Ha = \g\J^ — AJ3 can be diagonalized in a direct way. 
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Josephson and Rabi regimes. In these cases 1/J^<<r<<l and 1 << 
r, respectively. Based on the above formulas, one finds (J3) ~ JrA/2Vo 
(Josephson case) and (J3) ~ JA/2Vo (Rabi case) giving the mean dimen- 
sionless AM per boson 



1 + 



tA 

2Vn 



A 

2Vo 



respectively. Owing to formulas ^1 and El in the Josephson case, the range 
of parameter A is [— 2J|5f|, < 2J|(7|]. For this regime, the further condi- 
tion [HI reduces to (2nri/VJ)i/2 + (A/2J|^|) < 1. In the Rabi case, con- 
dition [inion (J3) entails that A ranges in [—2JT\g\, 2JT\g\] {2JT\g\ = 2Vo) 
which is, in principle, much larger than the range allowed in the Joseph- 
son case. Considering once more condition E[ this gives in the Rabi case 
(2nr^/^/J)^/^ + (A/2V0) < 1. On easily checks that weakly excited states 
\En) satisfy the conditions on the restricted range of A provided n « J, 
and |A| << 2J\g\, |A| << 2Vo in the Josephson case and in the Rabi case, 
respectively. In both cases the latter inequalities represent condition in 
its stronger version. 



3 The Inonii - Wigner contraction in the re- 
pulsive case 

The classical study of repulsive Hamiltonian Hr = —{\g\ J| -|- 2VoJi + AJ3), 
discussed in|Xl shows that, with r = Vo/J\g\ > 1 (Rabi regime), the energy 
minimum is such that Ji = J, J2 = J3 = 0. As shown by equation EHl 
a generic state near the minumun is such that Ji ~ J, IJ2I, IJ3I ^ J. 
In the Fock/ Josephson regimes, where r = VQ/J\g\ < 1, Hamiltonian Hr 
displays two minimum-energy states (see equation I^Uj) entailing low-energy 
configurations characterized by J3 ~ ±J, | J2I, |Ji| <^ J. 



3.1 Repulsive regime with r > 1 

In the Rabi regime (r > 1), the CPA valid for the attractive model can be 
implemented again. Then assuming hi, h2, hs as in formulas El the result of 
the contraction gives Ji — >• J — n, J2 — —^/Jp, and J3 \^q, which reduce 
Hr to a quadratic form. By defining Q = q — c, with c = \/JA/ {2VoW'^), the 
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final form of Ha is found to be 



Vo 



JA2 



(17) 



Since the eigenvalues of + l^^g^ are 
Hr is 



Vo 



2W{n + l/2) -2J 



2W{n + 1/2), the spectrum of 
JA2 



(18) 



As in the attractive case, the eigenfunctions $„((5) of Hamiltonian IT7I allow 
one to determine components Xm through the formula Xm{En) = $n(<5)- 
The energy eigenstates turn out to be 



\En) = T.mXm{.En)\J]m) , Xm{E„ 



WQ)e 



wet 



(19) 



This description is valid if the conditions on the 



with Q = mjyj — c. 
gaussian deviation and the Hermite-polynomyal zeros < J2/W and 

IcI < \fj - 



^2n/W , respectively, which can be rewritten as 1 < 2 Jr/ (r — 1) 
and |A|r/[2Vo(r-l)] < l-{2nri/V[J(r-l)i/2]}i/2^ are satisfied. For r > 1 , 
the first condition is fulfilled, while the second one gives A/2Vo < 1 — \j2n/ J . 
The latter is satisfied if A/2Vo < 1. Weakly excited states \En) with n > 
can be also considered provided J ^ 2n. Under such conditions, the mean 
dimensionless AM per boson is a linear function of A 



;i+(^3)/j) 



1 + 



Ar 



2Vo(r - 1) 



( J3) = cTJ = JAt/2Vo{t - 1) 



giving {Iz) ~ (1 + A/2Vo)/2 for r >> 1. Notice that (J3) coincides with 
formula IH^ for the minimum of the classical repulsive model and that, in the 
Rabi regime, (Z^) has the same form both for attractive bosons {g < 0) and 
for repulsive bosons ((7 > 0). 



3.2 Repulsive case with r < 1 

In this case, the classical ground-state configuration corresponds to two min- 
ima. The contraction scheme can be implemented in two ways by assuming 
/i2 = xJ2, hi = xJi, and /is = J3 =F I/x"^ which entails [hi,h2] = iix^h^, 
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[/i2, /is] = ihi, and [/13, hi] = ihi. Notice that /^s = J3 =F I/x"^ allows to de- 
scribe the two classical minima by further selecting a suitable definition for 
hs- The result of the contractions demonstrates the two possible choices 

hs = -n, h2 = J2/ "/j ^ p, hi = Ji/ VJ q, (20) 

and 

hs = +n, h2 = J2/VJ -p, hi = Ji/Vl q, (21) 

that are naturally associated to the Ja-positive and Js-negative minimum, 
respectively. The repulsive Hamiltonian Hr = —{\g\ Jl + 2VoJi + AJ3) thus 
can be cast in the two (local) forms 

= -\g\ J2 _ 2Jn + 2r + —(J - n) , (22) 

^ 1^1 ^ 

where s = ± recalls the presence of two minima. Notice that Hr could 
be diagonalized by means of the procedure used in the attractive case, pro- 
vided one adopts the rotated basis {\m)i = exp(— z7rJ2/2)|m)} of Ji and 
regards the q eigenvalues m/ ^/j as a continuous index. Unfortunately, while 
the evaluation of the energy eigenvalues is very easy in the "rotated" Ji 
basis {\m)i = exp(— i7rJ2/2)|m), |m| < J}, the eigenstates must be counter- 
rotated to recover the Js-basis representation that we have adopted in the 
other cases/regimes. This is a difficult problem in that recovering the eigen- 
states description in the J3 basis requires that the transformation matrix 
element (m'| exp(— z7rJ2/2)|m) is calculated explicitly and is formulated in 
the limit where m/y/j is a continuous index. 

To skip this problem, we observe that, owing to formulas 12 ( )l and HTl derived 
by the contraction procedure, + J2 + Jf = J{J + 1) — can be rewritten 
as J| ~ 2Jn — J2 while J3 = ±(J — n). We thus obtain the linearized 
expression J| ^ — ± 2JJ3. Hamiltonian 1^ reduces to Hr = —\g\[—J^ ± 
2JJ3 + 2tJJi + (A/|5f|) J3], whose digonalization is rather simple owing to 
the linear dependence on su(2) generators. Rewriting the latter as H^ = 
\g\ =F (2J ± 5) J3 — 2tJJi where 6 = A/\g\, the unitary transformations 
U± = exp(=FiJ20±) entail 

H^ = \g\[j'TR±U±J^U+], (23) 
with R± = (2 J ± (5)2 -|- Ar'^.P. The action oiU± is given by 
U-J^Ul. = J3 cos 0_ — Jl sin 0_ , f/+ Jst/^ = J3 cos 0+ + Ji sin 0+ , (24) 
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where angles 0± are definded by tg0_ = 2r J/ (2 J — S), tg0+ = 2r J/ (2 J + 5). 
The energy spectrum is thus represented by the eigenstates and the eigen- 
values 



\E^) = U±\m), E^ = \g\[j^Tm^{2J±6y + 4T^J^\, (25) 

respectively. One should recall that, within the present approximation scheme, 
these eigenvalues are significant for \m\ ^ J. Moreover, we notice that 
[/± — > 1 for r ^ thus reproducing the correct spectrum of the uncoupled 
model. The eigenvalues corresponding to the energy minima are obtained by 
setting m = — J and m = + J for H~ and H^, respectively, and read 



:= = \g\ [J' - J^{2J ± 6)' + At' \. (26) 

The choice of the signs ±, and thus the recognition of the lowest-energy 
states, is related to the sign of 6. This is discussed below. The states asso- 
ciated with eigenvalues 1261 take the form of su(2) coherent states [22]. The 
standard su(2) picture of such states, also known as Bloch states, is given by 

|_j.j) = ,5..-r.-|_j) = |.^^dW> (27) 



with Cjs = y (2J)!/s!(2J — s)\, while the coherent-state labels z = \z\e^^ and 
C = are such that \z\ = tg|^|, z e C. Since the minimum-energy states 
have the form 

\E^) = e'^''"^^\±J), (28) 

where =piJ20± = =f(0±/2)(J+ — J-), the link with the coherent-state picture 
is almost immediate. Upon setting ^ = ^(f)±/2, the corresponding z reads 
z = =ptg((/)±/2) = =p2rJ/(2J± 6). In view of this, eigenstate l-E^) takes the 
new form 

2J 

=cos2^(0_/2)^Cj,tg^(0^/2)|s- J). (29) 

s=0 

If 5 < 0, state \E^j{6)) (we make explicit the dependence from 6 to illus- 
trate clearly the difference between the absolute minumum and the local 
minimum) corresponds to the lowest-energy state with eigenvalue Ej.j{6) = 
\g\[.P - J^{2J +\5\y + At'], since Em{5) < Eti{5) (see equation |2SI)- The 
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Figure 2: Both panels concern the repulsive case. Grey (dark) squares, 
diamonds, points and (piecewise-linear, continuous, dotted or dashed) curves 
are relevant to = 20 (A^ = 40). Left panel: r = 0.6 (Josephson regime), 
z/ = 0.8. Within the CMA, the ground-state components X^s, given by 
formula 1221 and described by squares, well approximate the X^s (edges of 
the grey/dark piecewise-linear curves) calculated numerically. Right panel: 
r = 1.6 (Rabi regime), z/ = 0.8. Points (diamonds) -joined by dashed/dotted 
lines to better distiguish different cases- describe ground-state X^s within 
the CMA (CSSA) referred to formula ITUl f formula IHHj) . Continuous piecewise- 
linear curves represent X^s obtained numerically. The CSSA is qualitatively 
better than the CMA approximation where curves are shifted on the right. 
Further comments are given in section El 
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remaining state \E^{S)) represents the local minimum found in the classi- 
cal dynamics. In the opposite case 6 > 0, the lowest energy state iden- 
tifies with \E^.j). This in fact corresponds to (see equation OBj) E~l^{5) = 
\g\[J^ - J^(2J + 5)2 + 4r2], which satisfies E^^{6) < Em{6) for 6 > 0. 
Notice that £'^j(— = E^{6). This feature is important because it con- 
firms the symmetry property e^^'"'^^Hr{5) e~*'^"'i = Hr{—S) of repulsive Hamil- 
tonian Hr{6) = —\g\{Ji + 2JrJi + 6J3) stating that the spectra of the 
cases 6 > and 6 < must coincide, the relevant Hamiltonians being 
related by a unitary transformation. Based on this fact, we find as well 
\Etj{6)) = e+'''-^^\E]^,j{-\5\)). By acting with e^'^^i on \E^{-\6\)) we get the 
expression 

\E+{6)) = e-^-^2</.-g+.-Ji I - J) = e*J-e-*J2*-| (39) 

[notice that 0_ = 0_(— where we have used the property of the Js-basis 
states e^'^^'"\m) = e^'^'^l —m). Upon observing that 0_(— = 0+(+|(5|) we 
conclude that the unitary transformation reproduces, up to a phase factor, 
the diagonalization-process formula = e-'-^^'^+ \ + J) in a consistent 

way. Therefore, the ground state of the case 5 > is obtained by calculating 
formula inni explicitly, which gives 

2J 

\E+) = cos'' (0+/2) ^ CjM' (0+/2) \ J-s), (31) 

where (/)+(+|5|) = 0_(— 15|). We notice that \E^) corresponds to a coherent 
state \+J,0 = e^"'+~^*'^~| +J) whose extremal state is |J) (instead of | — 
J)) where \v\ = tg\^\ with v = —tg{(j)_/2) reproduces ED As in the case 
6 < 0, the remaining state \E^j{6)) describes the quantum counterpart of 
the local minimum. The expectation value of J3 is easily carried out. By 
using equations one finds {{Jk)± = (-^mI'^/cI-^'m)' ^ — 1)2,3) {J3)± = 
(±J|(J3 cos0± =F sin0±) I ± J), (Ji)± = Jsin0±, (J2)± = 0, namely 

(^3)± = ^=ii^, W± = 7==, (32) 

where fi± = 2tJ / (2 J ± 5), which, expanded up to second order in r, appear 
to be consistent with the classical values ^2 of the minimum-energy config- 
urations. The choice + (— ) for the lowest-energy state, corresponding to 
5 > (5 < 0), entails 2J ±5 = 2J + \5\ in /i±. Thus (73)+ and {J^). simply 
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Figure 3: Repulsive case, r = 2.4 (Rabi regime), z/ = 0.8. Grey (dark) 
diamonds, points and piecewise-linear curves are relevant to = 20 (A^ = 
40). Diamonds (points) -joined by dashed/dotted lines to better distiguish 
different cases- describe ground-state X^'s given by formula IT^ fformula IHHjl 
within the CMA (CSSA). Piecewise-linear curves have the usual meaning. 
Both CSSA and CMA are satisfactory. Further comments are given in section 

E 

differ of a factor —1. In passing we notice that states \E^f{6)) with the same 
6, should satisfy the condition {E^.j\E^.j) = they corresponding to different 
eigenvalues. \E^{5)) obtained within the CPA can be shown to be almost 
orthogonal ^T]. Excited states labeled by m = ±(n — J) with n « J 
can be derived explicitly from formula 123 By expressing them as \E^) = 
f/±J^| ± J)/(n!Cj„), one obtains = (U^J^Ul T tJ2r\E^j) /(nlC.^) , 

with U±JiU± = cos(0±)Ji =Fsin(0±)J3, that can be used to calculate the ex- 
pectation values of operators Jk, = 1, 2, 3 . The condition under which the 
eigenvalue that corresponds to the local minimum represents the first excited 
state can be determined quite easily (e. g., for 6 < 0) from E^{S) < E~{S) 
with m = — J + 1. 

Within Fock and Josephson regimes (t < 1), the AM per boson is readily 
evaluated from formulaElgiving (4) = [1±{2J±S) / ^Ar^J"^ + {2J ± 5y]/2. 
If r << 1, due to 0± ~ 2r J/ {2 J ±6) and in view of equations l29l and 130} the 
ground state reduces to \E^) ~ [1 - 2J{(j)^/2f][\ T J) + T ^ ± 1)] 
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where — and + are related to the cases S < and S > 0, respectively. 
Thus in the Fock regime (r << it is natural to set 0± ~ 0. By 

neglecting also the first order corrections, the ground state is approximated 
by \Eg{6)) = 9{6)\J) + 9{—S)\ — J) which, inserted in formula 221 gives 
(4) = e{5) = {lT l)/2. This well matches the case r = where Eg{±\5\) = 
-l\g\J^±JA) with5 = A/|^|. 

4 The coherent-state semiclassical approxi- 
mation. 

An alternative way to approximate both the ground state and the correspond- 
ing energy is to find the quantum counterpart of a classical configuration in 
terms of coherent states. If the hamiltonian algebra of a given model is 
known together with the coherent state relevant to such an algebra, classi- 
cal variables can be put in a one-to-one correspondence with the complex 
labels parametrizing a coherent state [22j. This is the case for Hamiltonian 
EJandlHthat are written in terms of su(2) generators J3, J±. Coherent states 
| — J, ^) of algebra su(2) are defined by equation |2ZI The latter allows one to 
parametrize a coherent state by z since ^ = |^|e*^ is related to z = |z|e*^ by 
kl = tg|^|. For a generic \z) the expectation values (J^) = {z\Jk\z), k = ±,3, 
given by 

{Js) = A\z\'-l)/i\z\' + l), {.U) = 2Jz*/i\z\' + l), (33) 

with (J_) = allow one to determine z when (J^) are known. Notice 

that (Ji) = ((J+) + (J-))/2 and (J2) = ((J+) - {J-))/2i. Therefore clas- 
sical configurations characterized by known values of Ji, J2 and J3 can be 
associated with a specific z by identifying each classical Jk with (J^) and 
observing that, owing to equations EH the phase 9 of z coincides with the 
phase of J+ = Ji + iJ2 while |zp = {J + J-i)/{J — J3). Recalling that 
this assumption becomes exact in the semiclassical limit J — > cxd, we name 
the map Ji, J2, J3 z coherent-state semiclassical approximation (CSSA). 
Determining J^'s that characterize the classical energy minimum thus pro- 
vide the ground-state approximation \Em) — \z) where \z) is determined 
by the previous semiclassical map. The correponding energy is obtained by 
i?- = {z\H\z). 
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5 Conclusions 



We have discussed the effectiveness of the CPA based on the Inonii-Wigner 
transformation by comparing the ground state (GS) obtained in the vari- 
ous regimes of both the repulsive and the attractive models with the ex- 
act lowest-energy eigenstate determined numerically. In the attractive case 
{g < 0), both for r < 1 and for r > 1, and in the repulsive case {g > 0) for 
r > 1 the CPA leads to approximate X^'s of weakly excited states through 
the eigenfunctions of equivalent harmonic-oscillator problems represented by 
formulas^! and El respectively. Due to the presence of two classical minima 
in model EJ the repulsive case with r < 1 requires that a different diagonal- 
ization scheme is developed after implementing the CPA on Hamiltonian 
lU This involves weakly excited states represented in terms su(2) coher- 
ent states EHl and EU In the attractive case, figures [T] show that the exact 
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0.5 1 1.5 2 2.5 3 

X 

Figure 4: Attractive case. Comparison of the exact ground-state (GS) energy 
with the ground-state energy El within CPA. Different colors in the tu plane 
are related to different value of indicator a. White regions are characterized 
by an excellent agreement of the exact and the approximated GS energies. 
See section El for details. 

components (calculated numerically) are almost indistinguishable from com- 
ponents Xm's obtained within the CPA and described by formula El The 
cases = 20 and = 40 that correspond to r = 0.02 > l/J^ = g.Ol and 
T = 0.02 > 1/J^ = 0.0025, respectively, describe the approach from above to 
the lower bound of Josephson regime. In the repulsive case, figures |21 allow 
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Figure 5: Repulsive case. Comparison of the exact ground-state (GS) energy 
with approximate GS energies. Different colors in the tv plane are related 
to different value of indicator a. White regions are characterized by an 
excellent agreement of the exact and the approximated GS energies. Details 
are discussed in section |S1 Left panel: a for the GS energv ITHl ir > 1). 
Right panel: a for the GS energy within the CSSA. (expectation value of 
Hr for the GS relevant to formula IH^ . 



one to compare the exact components (calculated numerically) with compo- 
nents Xm's obtained within the CPA and described by formula 1221 for r = 0.6 
(Josephson regime), u = 0.8 and N = 20, 40, and by formula IT^ for r = 1.6 
(Rabi regime), u = 0.8 and = 20, 40. While in the first case formula EH 
representing a su(2) coherent state, provides a satisfactory approximation, 
in the second case formula ^1 exhibits a shift on the right of highest weight 
components X^'s that, in addition, are smaller than the exact ones. In figure 
121 (right panel) X^s evaluated within the CSSA better match the exact ones 
both qualitatively and quantitatively. When r is increased (see figure El), the 
CPA approximation (CSSA) is satisfactory even if it tends to underestimate 
(overestimate) exact X^s. Figures 01 and El illustrate, through the parameter 
a = {E%i — E'J^^)/AE, the deviation of the GS energies obtained within the 
CPA or the CSSA from the GS energy calculated numerically. Energies El^, 
E'^, and AE are the exact GS energy, the approximated GS energy and the 
energy range defined as AE = E^^^ — Elj, respectively. E^^^ is the exact 
maximum energy. White, light grey, and dark grey colors identify the regions 
in the r u plane where a < 0.001, 0.001 < a < 0.01 and 0.01 < a < 0.1, re- 
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spectively. In figure El describing the attractive case, E'J^ is given by formula 
El E'^ well approximates the exact GS energy in the large (white) region in 
the T u plane. The repulsive case is considered in figure El In the left panel, 
E'^ given by formula ITHl is shown to well approximate the exact GS energy 
in a rather restricted region in the tu plane. On the contrary, right panel 
shows that evaluating E'J^ based on ground-state ESI within the CSSA pro- 
vides the best approximation (a < 0.001) almost everywhere. Concluding, 
except for the repulsive Josephson regime, where the CPA is not satisfactory, 
both the CPA and the CSSA provide a satisfactory approximation. The CPA 
is particularly good in the attractive-boson case. Among the many applica- 
tions to bosonic-well systems currently studied, such approaches seem quite 
appropriate for studying the low-energy spectrum of the three-well boson 
systems where the complexity of the energy-level structure mirrors the dy- 
namical instabilities of the chaotic three- well classical dynamics |2S1- The 
study of similar aspects in the three- AM mode rotational fluid outlined in 
jS] is currently in progress. 



A Classical energy minima 

The classical version of the attractive model El displays a dynamics charac- 
terized by four (two) fixed points if 1 S> r (r ^ 1). This can be seen by 
considering the relevant motion equations 

ji = (A-2|(7|J3)^2, j3 = -2VoJ2, j2 = 2{\g\J^ + Vo)Js-AJ^, (34) 

equipped with the motion constant J"^ = Jl + + Jl, that entail the fixed- 
point equations J2 = 0, 2\g\J^Ji + 2V0J3 — AJi = 0, with the constraint 
■P = + Jf. Their exact solution involves a fourth-order equation in J3, 
except for A = when the possible solutions are either Ji = — Vo/|5'| = —Jr 
or J3 = 0. In the general case A 7^ 0, if 1 3> t and J\g\ > A > (namely, for 
A sufficiently small), the searched solutions are such that either J3 ~ ±J, 
J > I Ji I , or 

Ji-±J, J > J3 > . (35) 

This feature can be proved explicitly. Particularly, the second pair of solution 
is obtained by implementing the approximation Ji = _ j2 ^ sJ{l — 

J|/2 J^), s = ±1. Neglecting the third order terms in J3/ J, the second fixed- 
point equation becomes (A/2 J) J| + 2J|(7|(1 + sr) J3 — J A = 0, whose roots 
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are found to be J3 = 2Ja^^[-l ±^l + (yl/2] with a, = A/[J|^|(1 + sr)]. 
While the negative root must be discarded because it entails IJ3I > J, the 
positive root -this can be shown to describe both a minimum (s = +1) and 
a saddle point {s = —1)- can be approximated as 

JtA , , 

(36) 



" 2K)(l + sr)' 

if (5 = A/|(7| < J. When r > 1 (and thus for r ^ V) the choices s = —1, +1 
are related to a maximum and a minimum, respectively. Notice that the 
previous formula giving the J3 coordinate is well defined for the minimum 
[s = +1) also when r ^ 1. 

Let us consider now the (classical) repulsive model El The corresponding 
Hamiltonian equations read 

ji = (A + 2|(7|J3)^2, j3 = -2^J2, j2 = 2(1^0 -|^7|^i)^3-AJi, (37) 

and exhibit once more the motion constant J"^ = Jf + +J2 + J^- For A = 
and r > 1, the energy minimum is easily shown to correspond to Ji = J, 
J2 = J3 = 0. Thus a generic state near the minumun is such that 

Ji~J, IJ2I, |J3|<^- (38) 
If A 7^ 0, provided A/J|5f| is sufficiently small, this statement is certainly 



valid for 1 ^ r = Vq/ J\g\ (Rabi regime). In fact, by setting Ji = y — 
~ J(l — J|/2J^) and neglecting the third order terms in J3/J in the fixed- 
point equation = 2{Vo — \g\Ji)J3 — AJi, one finds (5/2 J) J| + 2 J(r — 1) J3 — 



5J = 0, whose roots are found to be J3 = 2Ja~^[— 1 ± yl + a^/2], with 
a = A/[J\g\{T — 1)]. Discarding the negative root which entails IJ3I > J, 
the positive root can be approximated as 

J3 - ^ - (39) 
'^'-2\g\{r-l)-2Vo{r-iy ^^^^ 

if A/J\g\ <^ r — 1. In the Rabi regime where 1 <C r ~ r — 1 such condition 
reduces to A ^ Vq. In the Fock/Josephson regimes, where r < 1, the two 
configurations Ji = tJ, J3 = ±Ja/1 — are found to minimize the energy 
if A = 0. This suggests that, even with A 7^ 0, low-energy states are such 
that 

J3^±J, IJ2I, |Ji|<J. (40) 
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To obtain the energy- minimum configurations, in addition to J2 = 0, we 
consider the second fixed-point equation under the approximation J3 = 



syJ^ — Ji ~ sJ{l — Jf/2J^) with s = ±1. Neglecting the third order 
terms in Ji/J, the latter entails = {Vo/J)Jf + {2\g\J - sA)Ji - 2VqJ , 
which supply, with s = +1, two minimum-energy configurations {5 = A/\g\) 



T 



J 



These reproduce correctly the formula of the case A = 0. 
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